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Counting Statistics in Nanoscale Junctions. 
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We present first-principles calculations for the third moment of the current in atomic-scale junc- 
tions. We calculate this quantity in terms of the effective single-particle wave-functions obtained 
self-consistently within the static density-functional theory. As an example, we investigate the rela- 
tions among the conductance, the second and third moments of the current for carbon-atom chains 
of various lengths bridging two metal electrodes. We find that the conductance, the second-order 
and the third-order Fano factors show odd-even oscillation with the number of carbon atoms with 
the third-order Fano factor positively correlated to the conductance. 



Nanoscale electronics has generated a tremendous 
wave of scientific interest in the past decade due to 
prospects of device-size reduction offered by atomic-level 
control of certain physical properties [l|. In addition, 
it has spurred great interest in the fundamental under- 
standing of quantum transport [21 . One of these funda- 
mental questions relates to the moments of the current. 
For instance, the second moment - shot noise - defines 
the quantum fluctuations of the current at zero tem- 
perature due to the quantization of charge. Shot noise 
reaches the classical limit 2e/, where e is the electron 
charge and / is the average current [3|, when electrons 
in a conductor drift in a completely uncorrelated way as 
described by a Poissonian distribution of current events. 
On the other hand, in the presence of a junction with 
a narrow constriction, electrons, shot noise can be ex- 
pressed as 5*2 oc X^n^" (1 ~ -^n) ^^ tcrms of the trans- 
mission probabilities of each eigen-channel T„ [J]. It is 
a powerful tool for the exploration of quantum statistics 
of non-equilibrium electrons [5l4l2l| and may provide a 
means to explore also local temperature effects in nano- 
structures [Il[. In fact, it has been recently employed to 
characterize the signature of molecules/atomic wires in 
junctions [a la. 1 13| . 

The higher moments of the current (although more dif- 
ficult to measure and calculate) provide deeper insight 
into the statistics of charge dynamics, and are there- 
fore more refined tools to characterize the signature of 
molecules in junctions. However, no studies have con- 
sidered higher moments of the current in truly atomic- 
scale systems. To address this issue, we have developed 
a theoretical approach that combined with static density 
functional theory (DFT), which allows us to compute cor- 
relations up to the third moment of the current in atomic 
junctions. 

We have then investigated the relation between the 
conductance, the second moment (or shot noise, 6*2), and 
the third moment of the current (we denote it with 53) for 
a prototypical nanojunction consisting of an atomic chain 
with different number of carbon atoms connecting two 
metal electrodes as shown in Fig. [Ifa). This is not just 
an academic example since carbon is a versatile element 
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FIG. 1: (color online) (a) Schematic of the four-carbon 
atomic junction, (b) The spacial distribution of partial charge 
density for electrons with energies near the Fermi levels shows 
TT-orbital characters at Vb = 0.01 V. Vertical black lines corre- 
spond to the edges of the jellium model and circles correspond 
to atomic position. 



capable of forming diverse structures including diamond, 
graphite, fullerenes, nanotubes and graphene. Recently, 
experimentalists have shown the possibility to form car- 
bon atomic chains from graphite using a transmission 
electron microscope IJ]. The carbon atom chains have 
regularly patterned electronic structures as a function of 
the number of carbon atoms [l5J. In this regard, they 
are among the few model systems in which theory and 
experiments can be reasonably compared. 

Let us then consider a system at steady state subject 
to a bias Vb = {Efr — EpL)/e, where EpB. and EpL are 
the right and left electrochemical potential, respectively. 
The system is described by the field operator 
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FIG. 2: (color online) The density of states in the contin- 
uum region for N-atom carbon chains sandwiched between 
two metal electrodes at Vb = 0.01 V. Zero of energy corre- 
sponds to the left Fermi level. 



where a ~ L or R; a^K (t) = cxp{—iujt)a^-^ , and 

'^EK ^^ ^^^ annihilation operators of electrons inci- 
dent from the left (right) reservoir, satisfying the anti- 
commutation relations 

{<,k,,4IkJ = SafiHEi - E2)5{K, - K2), (2) 
where /3 = L or R. The single-particle wave functions 
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K incident from the left (L) and right (R) electrodes 
are later computed self-consistently in the framework of 

DFT [i,[ii,[i7|. 

The current operator is defined as 
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At zero temperature, the average of current operator 
gives the first moment, 
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^K/^kVW: (4) 



< / >= — / dE dv 

where the following expectation values have been used: 

< <k,«1k. >= 5^p5 (El - E,) <5(Ki - K2)/|, (5) 

L( R) 

with f^ the Fermi-Dirac distribution function in the 
left (right) electrode. 

We note that our wave-functions as obtained in the 
framework of DFT calculations describe noninteracting 
electrons, where the effective single-particle wavefunc- 
tions have boundary conditions describing electrons that 
are partially transmitted and partially reflected. Unlike 
the formalism developed in, e.g., Ref. [4], which has in- 
out ordering, the current operator defined in Eq. (j4|) de- 
scribes steady-states current where the time ordering is 
not involved. Therefore, the current correlation functions 



are defined in the time-unordered way 2l| . 



S2iio) = 27Th / dih - i2)e'"(*^-*^) < Ai{h)Ai{t2) >, 



(6) 



, describe electrons with energy E and momentum and 
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where AI{t) = I{t) - (/); ^2(0;) and S3{uj,uj') are the state current can be defined as S2 
two- and three-current spectral functions, respectively. "JS — ■^sl'^ — U, w — U), 

The zero- frequency 2nd and 3rd moment of the steady- 
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where the expectation values in Eqs. ([6]) and ([7|) have orem [18| . 
been calculated using the Wick-Bloch-De Doniinicis the- 
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where Ai denotes either creation or annihilation opera- 
tors and 77 = —1(1) for Fermions (Bosons). We note that 
alternative definitions of Fourier transform, e.g., the inte- 
grals with respect to (ti — ^2) and (^2 — ^3) in Eq. (|6]) and 
d?]) may lead to different parametrization of frequencies a; 
and u)' in the spectral functions. However, in the case of 
the steady-state current where uj = w' = 0, the zero- 
frequency current correlations are independent of the 
choices of Fourier transform. We also note that Eq. (|8]) 



leads to the relations S'f^ = S^^ = ~S^^ - 
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where S^^ = ^2(2: -J> cxd,z -> 00), S^^ = 6*2(2; -^ 
— 00, z -^ —00), 82^ = '52(2; — > —oo,z — >■ 00), and 
5^^ = 52(2; -^ 00, z — !► —00), which are a consequence 
of current conservation. Similarly, Eq. (jH]) leads to 
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For a single-channel tunnel junction with transmis- 
sion probability T, the first, second, and third moment 
of current are given by J oc T, 5*2 oc T(l — T), and 
53 oc -2T2(1 - T), respectively. The result of the 
unordered third moment is consistent with the results 
of time-unordered three-current correlations derived by 
other groups [19|, l20|. Finally, we define the second- 
and third-order Fano factors (which is dimensionless) 
in the small bias regime for steady-state currents as 
F2 = S2/{2eI) and F3 = S'3/[(2e)2j], respectively. As 
a direct result, for the single-channel junction G oc T, 
F2 oc (1 - T), and F3 oc -2T(1 - T), respectively. 

As an example, we have investigated the counting 
statistics in linear atomic chains formed by four to ten 
carbon atoms (denoted as C4 to CIO) bridging between 
two metal electrodes modeled as electron jclhum (r^ — 2). 
The distances between two neighboring carbon atoms are 
2.5 a. u., and the end atoms of the chain are fixed at 



1.4 a. u. inside the positive background edge of the elec- 
tron jellium as shown in Fig.[ljb) \l^. As it was found in 
Ref. |15J when the length of the wire is increased by one 
carbon atom, two electrons are added to the 7r*-orbital 
as shown in Fig. [2j The odd-numbered chains have a 
higher conductance due to half-filled 7r*-orbital while the 
even-numbered chains have a lower conductance due to a 
full-filled 7r*-orbital at the Fermi levels, as shown in the 
upper panel of Fig. [S] In the middle and lower panels 
of Fig. [3] we show the influence of the number of car- 
bon atoms on the second-order F2 and the third-order 
Fano factor F3, respectively, in the linear response regime 
{Vb = 0.01 V). We observe that F2 and F3 both display 
odd-even oscillation with the number of carbon atoms. 

In order to better understand the relation among the 
moments of counting statistics in the carbon atom chains, 
we have investigated the differential conductance (defined 
as G = dl/dV), the differential second-order Fano fac- 
tor [defined as F2 = {l/2e)dS2/dI], and the differential 
third-order Fano factor [defined as F3 = {l/2e)'^dS3/dI] 
for the C4 and C5 wires in the non-linear response regime. 
We observe that the differential conductance of C4 chain 
increases as the applied bias increases, while the differ- 
ential conductance of C5 chain decreases as the applied 
bias increases, as shown in the top panels of Figure |4l( a) 
and (b), respectively. The increase (decrease) of dif- 
ferential conductance with the biases for the C4 (C5) 
chain is due to the half-filled (full-filled) n* orbital at the 
Fermi level, where more (less) states are included in the 
current-carrying energy window created by increasing bi- 
ases. The bottom panels of Fig. HJ^a) and (b) show that 
conductance and F3 are strongly positively correlated in- 
dicating that the dominant eigen-channels for counting 
statistics have transmission probabilities T > 0.5. 
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increases (decreases) as bias increases in even- (odd-) 
numbered carbon atom chains. We observe that F3 and 
differential conductance are significantly positively corre- 
lated, thus showing that third-order Fano factor provides 
more information than the second-order Fano factor re- 
garding the transmission probabilities of eigen- channels. 
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FIG. 3: (color online) Conductance [(black) square; top 
panel], the second-order Fano factor [(red) circle; middle 
panel], and the third-order Fano factor [(green) triangle; bot- 
tom panel] of the atomic wires as a function of the number of 
carbon atoms in the wire at Vb ~ 0.01 V. 




FIG. 4: (color online) (a) For four-carbon atom wire and 
(b) for five-carbon atom wire: the differential conductance 
[(black) square; top panel] , the second-order differential Fano 
factor [(red) circle; bottom panel], and the three-order differ- 
ential Fano factor [(blue) triangle; bottom panel] vs. bias. 



In conclusion, we have calculated the third moment of 
the current within DFT that allows the study of count- 
ing statistics at the atomic level. As an example, we have 
investigated the relation among conductance, second and 
third moments of the current for carbon chains of differ- 
ent length. In the linear response regime, conductance, 
second and third moments show odd-even oscillations 
with the number of carbon atoms, which is mainly due 
to the orderly patterned electronic structure of carbon- 
atom chains. In the nonlinear regime, the conductance 
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